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Objective:  To  test  the  hypothesis  that  components  of  heart 
period  variability  track  autonomic  function  during  simulated  hem¬ 
orrhage  in  humans. 

Design:  Prospective  experimental  laboratory  intervention. 

Setting:  Human  physiology  laboratory. 

Subjects:  A  total  of  33  healthy,  nonsmoking,  volunteer  subjects 
(23  men,  ten  women). 

Interventions:  Progressive  lower  body  negative  pressure  was 
applied  in  5-min  stages  until  the  onset  of  impending  cardiovas¬ 
cular  collapse. 

Measurements  and  Main  Results:  The  electrocardiogram,  beat- 
by-beat  finger  arterial  pressure,  and  muscle  sympathetic  nerve 
activity  from  the  peroneal  nerve  were  recorded  continuously. 
Pulse  pressure  was  calculated  from  the  arterial  pressure  wave¬ 
form  and  used  as  an  estimate  of  relative  changes  of  central  blood 
volume.  Heart  period  variability  was  assessed  in  both  time  and 
frequency  domains.  Application  of  lower  body  negative  pressure 
caused  progressive  reductions  of  R-R  interval  and  pulse  pressure 


and  progressive  increases  of  muscle  sympathetic  nerve  activity. 
Arterial  pressures  changed  minimally  and  late.  R-R  interval  time 
domain  variability  measures  and  spectral  power  at  the  high 
frequency  (0.15-0.4  Hz)  decreased  progressively  with  lower  body 
negative  pressure  (p  <  .001).  Both  R-R  interval  high-frequency 
power  and  time  domain  variability  measures  correlated  inversely 
with  muscle  sympathetic  nerve  activity  and  directly  with  pulse 
pressure  (all  amalgamated  R2  >  .88,  all  p  <  .001). 

Conclusions:  Components  of  heart  period  variability  track  early 
compensatory  autonomic  and  hemodynamic  responses  to  pro¬ 
gressive  reduction  in  central  blood  volume.  Such  analyses,  inter¬ 
preted  in  conjunction  with  standard  vital  signs,  may  contribute  to 
earlier  assessments  of  the  magnitude  of  blood  volume  loss  during 
hemorrhage.  (Crit  Care  Med  2008;  36:1892-1899) 

Key  Words:  central  hypovolemia;  lower  body  negative  pressure; 
stroke  volume;  pulse  pressure;  blood  pressure;  remote  medical 
monitoring 


Hemorrhage  is  one  of  the 
leading  causes  of  death  for 
both  civilian  trauma  pa¬ 
tients  and  battlefield  casual¬ 
ties  (1).  Despite  this,  first  responders 
have  no  tools  to  assess  volume  status 
before  definitive  medical  intervention. 
Standard  vital  signs  such  as  arterial  pres¬ 
sure,  heart  rate,  and  arterial  oxygen- 
hemoglobin  saturation  remain  stable  (if 
somewhat  high  in  the  case  of  heart  rate) 
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during  hemorrhage,  until  the  onset  of 
cardiovascular  decompensation,  and 
therefore  provide  little  information  on 
early  changes  of  volume  status  and  sub¬ 
sequent  identification  of  the  need  for 
prompt  intervention  in  severely  injured 
patients  (2). 

Importantly,  compensation  to  severe 
hemorrhage  is  driven  by  the  autonomic 
nervous  system.  Loss  of  blood  volume 
triggers  withdrawal  of  parasympathetic 
and  activation  of  sympathetic  neural  ac¬ 
tivity  to  defend  against  decreases  of  arte¬ 
rial  pressure  (2).  Consequently,  arterial 
pressure  may  fail  to  provide  an  early  in¬ 
dication  of  hemorrhage  severity  because 
it  is  maintained  fairly  constant  during  the 
early  stages  of  mild  to  moderate  hemor¬ 
rhage  through  autonomic  neural  mecha¬ 
nisms  that  increase  peripheral  vascular 
resistance  and  heart  rate.  By  the  time 
arterial  pressures  fall  to  dangerous  levels 
(which  can  occur  quickly),  it  may  already 
be  too  late  to  apply  effective  lifesaving 
interventions.  Instead  of  monitoring  reg¬ 
ulated  variables  (i.e.,  outcomes),  such  as 
arterial  pressure,  medical  personnel  need 
the  capability  to  monitor  responses  of 


controlling  mechanisms — in  this  case, 
autonomic  neural  function. 

Frequency  domain  analysis  of  heart 
period  variability  has  been  used  as  a  tech¬ 
nique  to  assess  autonomic  function  in 
patients  monitored  in  an  intensive  care 
unit  (ICU),  and  differences  in  frequency 
characteristics  were  found  to  be  associ¬ 
ated  with  patient  outcome  (3,  4).  Norris 
et  al.  (5)  have  demonstrated  the  ability  of 
time  domain  measures  of  heart  rate  vari¬ 
ability  to  accurately  predict  death  from 
traumatic  injury  up  to  12  hrs  before  it 
occurs.  Grogan  et  al.  (6)  have  suggested, 
based  on  thousands  of  individual  ICU  pa¬ 
tient  records,  that  volatility,  a  derivation 
of  beat-by-beat  heart  rate  variability  in 
which  heart  rate  is  sampled  every  1  to  4 
secs,  shows  promise  as  a  new  vital  sign 
for  ICU  patient  monitoring.  We  have 
shown  previously,  through  frequency- 
domain  analysis  of  heart  period  (R-R  in¬ 
terval)  variability  in  prehospital  trauma 
patients,  that  the  high-frequency-to- 
low-frequency  ratio  derived  from  a  Fou¬ 
rier  transform  (HF/LF)  is  higher  in  pa¬ 
tients  who  eventually  die  compared  with 
patients  who  survive  traumatic  injury  (7, 
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8).  These  differences  are  detectable  at  a 
point  in  time  when  standard  vital  signs 
such  as  mean  arterial  pressure,  heart 
rate,  and  arterial  oxygen- hemoglobin 
saturation  are  not  different  between  pa¬ 
tient  groups. 

If  some  component  of  heart  period 
variability  does  in  fact  track  autonomic 
compensation  to  reduced  central  blood 
volume  (hemorrhage),  it  would  be  an  at¬ 
tractive  candidate  for  use  in  noninvasive 
monitoring  systems  for  trauma  patients 
and  combat  casualties.  However,  scien¬ 
tists  and  clinicians  are  split  as  to  whether 
heart  period  variability  does  (9,  10)  or 
does  not  (11,  12)  reflect  modulation  by 
the  autonomic  nervous  system.  Pagani  et 
al.  (13)  introduced  the  concept  that  the 
ratio  of  LF/HF  R-R  interval  spectral 
power  calculated  from  frequency  domain 
analysis  provides  an  indication  of  sympa¬ 
thetic  vs.  parasympathetic  cardiac  con¬ 
trol.  This  contention  derives  from  obser¬ 
vations  that  R-R  interval  low-frequency 
oscillations  (RRILF;  0.04-0.15  Hz)  are 
modulated,  importantly  (but  not  exclu¬ 
sively),  by  sympathetic  neural  traffic  (10, 
13).  However,  both  Kingwell  et  al.  (14) 
and  Saul  et  al.  (15)  failed  to  demonstrate 
significant  correlations  between  directly 
measured  muscle  sympathetic  nerve  ac¬ 
tivity  (MSNA)  and  RRILF  in  research  vol¬ 
unteers  at  rest.  Furthermore,  induced 
sympathetic  activation  with  nitroprusside 
infusion  sufficient  to  decrease  arterial 
pressures  by  about  15  mm  Hg  revealed 
nonsignificant  correlations  between 
MSNA  and  RRILF  (15).  Although  this  ob¬ 
servation  suggests  a  lack  of  association 
between  these  variables,  pharmacologic 
reduction  of  arterial  pressure  does  not 
represent  the  same  physiologic  condition 
as  the  progressive  reduction  of  central 
blood  volume  that  accompanies  severe 
hemorrhage.  However,  severe  hemor¬ 
rhage  in  humans  cannot  be  studied  in  a 
controlled  laboratory  environment. 

We  have  proposed  the  use  of  lower 
body  negative  pressure  (LBNP)  as  a  re¬ 
search  tool  to  simulate  the  progression 
from  mild  to  severe  hemorrhage  in  hu¬ 
mans  during  the  preshock  state  (16).  It 
has  been  known  for  some  time  that  pro¬ 
gressive  LBNP  reduces  heart  rate  vari¬ 
ability  (17,  18)  and  increases  MSNA  (19). 
However,  a  robust  assessment  of  the  di¬ 
rect  relationship  between  heart  rate  vari¬ 
ability  and  MSNA  during  progressively  se¬ 
vere  reductions  in  central  blood  volume 
(LBNP  levels  of  >45  mm  Hg)  has  been 
limited  because  recordings  of  sympa¬ 
thetic  traffic  (recorded  most  commonly 


from  tibial  or  peroneal  nerves  at  the 
knee)  are  often  lost  due  to  interference  of 
the  microelectrode  by  increases  of  nega¬ 
tive  chamber  pressure.  Negative  pres¬ 
sures  on  the  order  of  —60  mm  Hg  and 
above  induce  physiologic  responses  sim¬ 
ilar  to  blood  loss  of  >1000  mL  (16). 
Therefore,  the  issue  of  whether  changes 
of  heart  period  variability  reflect  in¬ 
creases  of  MSNA  during  severe  central 
hypovolemia  remains  to  be  resolved  be¬ 
fore  heart  period  variability  can  be  con¬ 
sidered  as  a  diagnostic  tool  for  assessing 
severity  of  hemorrhage. 

To  address  this  issue,  we  measured 
MSNA  from  peroneal  nerves  directly,  us¬ 
ing  the  microneurography  technique,  in 
a  cohort  of  healthy  human  volunteers 
during  simulated  hemorrhage  with 
LBNP.  The  purpose  of  this  study  was  to 
test  the  hypothesis  that  components  of 
heart  period  variability  track  autonomic 
function  during  progressive  central  hypo¬ 
volemia  of  a  magnitude  estimated  to  be 
similar  to  that  of  severe  hemorrhage. 

MATERIALS  AND  METHODS 

Subjects.  A  total  of  33  nonsmoking  sub¬ 
jects  (23  men,  ten  women)  volunteered  to 
participate  in  this  study  (age,  31  ±  2  yrs; 
height,  174  ±  2  cm;  weight,  80  ±  3  kg). 
Before  inclusion,  all  subjects  underwent  a 
medical  history  and  physical  examination  by  a 
physician  to  ensure  they  had  no  previous  or 
current  medical  conditions  that  might  pre¬ 
clude  their  participation.  To  ensure  female 
subjects  were  not  pregnant,  all  potential  fe¬ 
male  subjects  underwent  a  urine  test  within 
24  hrs  of  experimentation.  All  subjects  main¬ 
tained  their  normal  sleep  patterns,  refrained 
from  exercise,  and  abstained  from  caffeine  and 
other  autonomic  stimulants,  including  pre¬ 
scription  or  nonprescription  drugs,  for  >24 
hrs  before  the  study.  Subjects  received  a  verbal 
briefing  and  written  descriptions  of  all  proce¬ 
dures  and  risks  associated  with  the  study  and 
were  made  familiar  with  the  laboratory,  the 
protocol,  and  procedures.  Subjects  were  en¬ 
couraged  to  ask  questions  of  the  investigators, 
and  then  they  signed  an  informed  consent 
form  that  had  been  approved  by  the  Institu¬ 
tional  Review  Board  for  the  protection  of  hu¬ 
man  subjects  in  research  from  Brooke  Army 
Medical  Center  and  the  U.S.  Army  Institute  of 
Surgical  Research,  Fort  Sam  Houston,  TX. 

Experimental  Protocol.  Subjects  were  in¬ 
strumented  with  a  standard  four-lead  electro¬ 
cardiograph  (ECG)  to  record  cardiac  electrical 
potentials,  and  a  finger  cuff  to  record  beat-by- 
beat  finger  arterial  pressure  (Finometer  Blood 
Pressure  Monitor,  TNO-TPD  Biomedical  In¬ 
strumentation,  Amsterdam,  The  Netherlands). 
MSNA  was  recorded  directly  from  peroneal 
nerves  with  the  microneurography  technique 


(20).  Central  hypovolemia  was  induced  by  ap¬ 
plication  of  LBNP  to  simulate,  as  closely  as 
possible  in  healthy  human  volunteers,  the  he¬ 
modynamic  challenges  associated  with  severe 
hemorrhage  (16).  Subjects  were  positioned 
supine  within  an  airtight  chamber  that  was 
sealed  at  the  level  of  the  iliac  crest  by  way  of  a 
neoprene  skirt.  Because  injured  patients  do 
not  breathe  at  a  fixed  rate,  we  did  not  attempt 
to  control  breathing  frequency,  and  subjects 
were  allowed  to  breathe  spontaneously.  Each 
subject  underwent  exposure  to  progressive 
LBNP  until  the  point  of  impending  cardiovas¬ 
cular  collapse.  The  LBNP  protocol  consisted  of 
a  5-min  control  period  (baseline)  followed  by  5 
mins  of  chamber  decompression  at  —15,  —30, 
—45,  and  —60  mm  Hg,  and  then  additional 
increments  of  —10  mm  Hg  every  5  mins  until 
the  onset  of  cardiovascular  collapse  followed 
by  a  10-min  recovery  period.  Cardiovascular 
collapse  was  defined  by  one  or  a  combination 
of  the  following  criteria:  a)  onset  of  presynco- 
pal  symptoms  such  as  gray  out,  b)  a  precipi¬ 
tous  fall  in  systolic  pressure  of  >15  mm  Hg  or 
a  sudden  bradycardia,  c)  progressive  diminu¬ 
tion  of  systolic  pressure  to  <70  mm  Hg,  and 
d)  voluntary  subject  termination  due  to  dis¬ 
comfort  from  symptoms  such  as  sweating, 
nausea,  or  dizziness. 

Data  Analysis.  Data  were  sampled  at  500 
Hz  and  digitized  to  computer  (WINDAQ, 
Dataq  Instruments,  Akron,  OH).  The  last  3 
mins  of  data  from  each  LBNP  stage  and  the  3 
mins  just  before  cardiovascular  collapse  were 
extracted  for  analysis.  Individual  R-waves  gen¬ 
erated  from  the  ECG  were  marked  at  their 
occurrence  in  time  and  used  for  subsequent 
identification  of  systolic  and  diastolic  pres¬ 
sures  generated  from  the  Finometer.  Pulse 
pressures  were  calculated  from  the  difference 
between  systolic  and  diastolic  pressure  and 
used  as  an  estimate  of  changes  in  stroke  vol¬ 
ume  (21).  All  analyses  were  performed  with 
commercially  available  data  analysis  software 
(WinCPRS,  Absolute  Aliens,  Turku,  Finland). 

For  representation  of  heart  period  variabil¬ 
ity  in  the  time  domain,  we  calculated  R-R 
interval  sd  (RRISD)  and  the  percentage  of  ad¬ 
jacent  R-R  intervals  that  varied  by  >50  msecs 
(pNN50).  For  representation  of  heart  period 
variability  in  the  frequency  domain,  R-R  inter¬ 
vals  were  replotted  using  linear  interpolation 
and  resampled  at  5  Hz.  Data  then  were  passed 
through  a  low-pass  impulse  response  filter 
with  a  cutoff  frequency  of  0.5  Hz.  Data  sets 
were  submitted  to  a  Fourier  transform  with  a 
Hanning  window.  The  magnitude  of  R-R  in¬ 
terval  oscillations  were  quantified  by  calculat¬ 
ing  the  power  spectral  density  for  the  signal 
(total;  0.04-0.4  Hz).  Signal  areas  were  sepa¬ 
rated  into  high-frequency  (RRIHF;  0.15-0.4 
Hz)  and  low-frequency  (RRILP;  0.04  -  0.15  Hz) 
bands.  To  account  for  potentially  large  inter¬ 
subject  variability  common  with  frequency 
analyses,  we  also  normalized  RRIHF  and  RRILF 
to  total  power  by  dividing  both  RRIHF  and 
RRIlf  by  the  total  power  (the  sum  of  RRIHF 
and  RRILF)  (22). 
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Figure  1.  Muscle  sympathetic  nerve  activity  ( MSNA )  from  a  peroneal  nerve  is  shown  for  one  subject  during  progressive  increases  of  lower  body  negative 
pressure. 


Statistical  Analysis.  A  one-way  analysis  of 
variance  for  repeated  measures  was  used  for 
comparison  of  dependent  variables  across 
treatments  (levels  of  LBNP).  A  two-way  anal¬ 
ysis  of  variance  for  repeated  measures  was 
used  for  comparison  of  dependent  variables 
across  gender  groups  (men  vs.  women)  and 
treatments  (levels  of  LBNP).  Duncan’s  mean 
separation  procedure  was  employed  for  levels 
of  LBNP  when  the  probability  of  LBNP  main 
effects  was  <.05.  Amalgamated  Pearson  cor¬ 
relation  coefficients  were  calculated  among 
variables  of  interest.  The  probabilities  of  ob¬ 
serving  chance  effects  due  to  our  experimental 
manipulations  are  presented  as  exact  p  values 
(23).  All  data  are  expressed  as  mean  ±  se 
unless  specified  otherwise. 

RESULTS 

Average  time  to  cardiovascular  col¬ 
lapse  was  1816  ±  60  secs,  revealing  that 
a  majority  of  subjects  experienced  symp¬ 
toms  late  into  the  -60  mm  Hg  stage  to 
early  into  the  -70  mm  Hg  (n  =  33).  No 
subject  reported  feeling  anxiety  or  dis¬ 
comfort  associated  with  LBNP. 

Application  of  LBNP  often  causes  the 
microneurography  electrode  tip  to  shift 
and  move  beyond  the  recording  region 
within  the  peroneal  nerve  fascicle.  Such 
minute  changes  in  electrode  position  de¬ 
grade  the  signal  in  some  subjects  to  the 


point  where  spontaneous  bursts  of  nerve 
traffic  are  no  longer  observable.  We  were 
able  to  maintain  good  MSNA  recordings 
at  lower  levels  of  negative  pressure,  but 
signal  quality  tended  to  degrade  as  cham¬ 
ber  pressure  increased.  For  this  reason  we 
did  not  quantify  total  MSNA  (burst  area 
multiplied  by  the  number  of  bursts)  but, 
rather,  quantified  MSNA  as  burst  rate 
(bursts  per  100  heart  beats).  Although  our 
protocol  took  each  subject  to  cardiovascu¬ 
lar  collapse,  we  only  included  data  from 
subjects  with  good  MSNA  recordings  for 
correlations  with  heart  period  variability 
metrics.  MSNA  was  maintained  throughout 
cardiovascular  collapse  in  12  subjects 
(seven  men  and  five  women),  and  there¬ 
fore,  the  unequal  sample  sizes  we  report 
for  each  LBNP  level  represent  primarily 
loss  of  MSNA  before  cardiovascular  col¬ 
lapse.  Of  the  12  subjects  in  whom  nerve 
recordings  were  maintained  until  cardio¬ 
vascular  collapse,  one  experienced  symp¬ 
toms  during  —60  mm  Hg,  six  during 
-70  mm  Hg,  and  five  during  -80  mm 
Hg.  An  example  of  an  MSNA  recording 
maintained  at  high  negative  pressures  is 
shown  for  one  subject  in  Figure  1.  MSNA 
increased  proportionately  with  increasing 
levels  of  LBNP  in  both  men  (F  =  8.391, 
p  <  .001)  and  women  (F  =  3.388,  p  = 


.011).  The  primary  physiologic  responses 
and  the  number  of  subjects  studied  at 
each  LBNP  stage  are  shown  in  Table  1. 

RRIlf  and  RRIHf  correlated  poorly 
with  MSNA  at  baseline  before  application 
of  LBNP  (r  =  .12  for  RRILF  and  MSNA; 
r  =  .10  for  RRIhf  and  MSNA;  n  =  33). 
Application  of  LBNP  caused  progressive  de¬ 
creases  of  R-R  interval  and  pulse  pressure 
but  changed  systolic  and  diastolic  pressures 
minimally.  Systolic  pressure  began  to  de¬ 
crease  at  —45  mm  Hg,  with  little  change  in 
either  diastolic  or  mean  pressure.  Respira¬ 
tory  rates  in  all  33  subjects  at  baseline 
(15  ±  1  breaths/min)  did  not  differ  statis¬ 
tically  (F  =  .532,  p  =  .600)  from  those  at  the 
maximal  level  of  LBNP  (15  ±  1  breaths/min) 
and  were  maintained  at  all  times  within  the 
high-frequency  range  (.15  to  .4  Hz). 

Application  of  LBNP  caused  progres¬ 
sive  reductions  of  RRISD,  pNN50,  and 
RRIhf  but  did  not  affect  RRILF  (p  =  .15). 
pNN50  and  RRIHF  were  reduced  at  -30 
mm  Hg  (both  p  <  .05),  whereas  RRISD 
was  reduced  at  -45  mm  Hg  (p  =  .04)  as 
compared  with  baseline.  These  patterns 
were  not  influenced  by  gender.  Progres¬ 
sive  reductions  in  RRIHF  were  observed  in 
both  men  (F  =  4.76,  p  <  .001)  and 
women  (F  =  3.99,  p  <  .001),  but  RRILF 
was  unaltered  across  LBNP  in  both 
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Table  1.  Hemodynamic  and  sympathetic  neural  responses  to  LBNP 


Variable 

Baseline 
n  =  33 

LBNP  15 
n  =  33 

LBNP  30 
n  =  24 

LBNP  45 
n  =  19 

LBNP  60 
n  =  14 

LBNP  70 

n  =  9 

LBNP  80 
n  =  5 

Recovery 
n  =  10 

RRI,  ms 

975  ±  23a 

951  ±  22b 

889  ±  28c 

804  ±  38d 

702  ±  47e 

633  ±  53f 

605  ±  55f 

1020  ±  58a,  b 

HR,  bpm 

63  ±  1.4a,  b 

64  ±  1.4a,  b 

70  ±  2.6a,  c 

77  ±  3.3c 

89  ±  5.4d 

99  ±  7.8e 

103  ±  8.5e 

61  ±  3.5b 

SAP,  mm  Hg 

131  ±  1.7a,  b 

132  ±  1.9a,  b 

127  ±  2.1b,  c 

122  ±  2.7c,  d 

119  ±  4.2d 

112  ±  8.1e 

116  ±  5. Id,  e 

136  ±  5.1a 

DAP,  mm  Hg 

77  ±  1.0a 

78  ±  1.1a,  b 

78  ±  1.1a,  b 

78  ±  1.2a,  b 

79  ±  1.8a,  b,  c 

77  ±  3.8a 

83  ±  3.0c,  d 

81  ±  2.8a,  b,  c 

MAP,  mmHg 

98  ±  1.2a 

98  ±  1.1a 

96  ±  1.1a 

95  ±  1.4a,  b 

94  ±  2.5b,  c 

89  ±  5.1d 

94  ±  3.2e 

102  ±  3.5a 

PP,  mm  Hg 

54  ±  1.4a 

54  ±  1.5a 

49  ±  1.9b 

44  ±  2.4c 

40  ±  3.2d 

34  ±  5.1e 

33  ±  4.6e 

53  ±  2.3a 

MSNA,  b/100 

21  ±  2.0a 

27  ±  2.3a,  b 

32  ±  2.6b,  c 

37  ±  3.5c,  d 

42  ±  3. Id 

41  ±  3.8d 

44  ±  5.7d 

26  ±  5.0a,  b 

Values  are  means  ±  se;  RRI,  R-R  interval;  HR,  heart  rate;  SAP,  systolic  arterial  pressure;  DAP,  diastolic  arterial  pressure;  MAP,  mean  arterial  pressure; 
PP,  pulse  pressure;  MSNA,  muscle  sympathetic  nerve  activity  expressed  as  the  number  of  bursts  per  100  heart  beats;  LBNP,  lower  body  negative  pressure. 
Differences  between  means  within  a  row  that  share  the  same  letter  are  not  distinguishable  statistically  (p  >  .05). 
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Figure  2.  Group  averages  are  shown  ±  SE  for  R-R  interval  SD  (RRISD),  the  percentage  of  adjacent  R-R  intervals  that  varied  by  >50  msecs  ipNN50),  and 
R-R  interval  oscillations  derived  from  a  Fourier  transform  at  the  low  ( RRILF )  and  high  (RRIHF)  frequencies.  Data  identified  with  brackets  show  differences 
from  baseline.  LBNP,  lower  body  negative  pressure. 


groups  (F  <  1.232,  p  >  .324).  There  were 
no  statistically  distinguishable  differ¬ 
ences  in  responses  of  heart  period  vari¬ 
ability  and  MSNA  between  men  and 
women  (F  <  1.134,  p  >  .442).  Changes  of 
time  and  frequency  domain  measures  of 
heart  period  variability  with  LBNP  for  the 
entire  group  of  33  subjects  are  shown  in 
Figure  2. 


After  normalization  of  R-R  interval  os¬ 
cillations  to  total  power,  RRILF  increased 
and  RRIHf  decreased  progressively  with 
LBNP  and  were  significantly  different 
from  baseline  at  -45  mm  Hg  ip  <  .001). 

Figure  3  shows  associations  among 
time  and  frequency  domain  measures  ex¬ 
pressed  in  absolute  units  vs.  MSNA.  In  both 
time  and  frequency  domains,  analyses  of 


components  of  heart  period  variability  dem¬ 
onstrated  high,  inverse  correlations  with 
MSNA  (R2  >  .95),  with  the  exception  of  a  low 
correlation  between  MSNA  and  RRILF  calcu¬ 
lated  in  absolute  units  (Z?2  =  .09). 

High  amalgamated  correlation  coef¬ 
ficients,  as  shown  in  Figure  3,  may 
sometimes  mask  a  high  rate  of  fairly  low 
individual  correlations.  To  address  this  issue, 


Crit  Care  Med  2008  Vol.  36,  No.  6 


1895 


35  -| 


30  - 


25  - 


z  15  - 

Q. 


10  - 


5  - 

0  — 
15 


R2  =  0.95 

I— O—l 

l— O—l 


li  Ti 

111  Ji 

20  25  30  35  40  45  50 


70  -| 


65  - 


60  - 
c n 

E  55  - 

Q 

W 

I  50- 
cc 

45  - 


40  - 


35  4- 
15 


T  R2  =  0.96 

MM 

HH 


20  25  30  35  40  45  50 


2400  i 
2000  - 

c,  1600  - 

m 

E 

lL  1200  - 

__j 

DC 

CL 

800  - 

400  - 

0  — 
15 


R2  =  0.09 


20  25  30  35  40  45  50 

MSNA,  b/100  hb 


cc 

CL 


1600 


1200  - 


800  - 


400  - 


0  4- 
15 


T  R2  =  0.98 

MM  " 

..  MM 


I — O — I 


20  25  30  35  40  45  50 

MSNA,  b/100  hb 


Figure  3.  Group  averages  for  time  and  frequency  domain  variables  ±  SE  are  plotted  against  muscle  sympathetic  nerve  activity  (MSNA),  with  amalgamated 
correlation  coefficients  shown  in  the  upper  right  corner  of  each  panel.  RRISD,  R-R  interval  SD;  RRILF,  R-R  interval  oscillations  derived  from  a  Fourier 
transform  at  the  low  frequency,  RRIhf,  R-R  interval  oscillations  derived  from  a  Fourier  transform  at  the  high  frequency;  pNN50,  the  percentage  of  adjacent 
R-R  intervals  that  varied  by  >50  msecs. 


Table  2.  Average  individual  correlations  among 
various  heart  rate  variability  metrics  and  MSNA 
during  LBNP  from  a  sample  of  subjects  from 
the  study  cohort  who  provided  data  up  to  -70 
mm  Hg  (n  =  9) 


Variables 

Mean  Pearson  r 

%  individual 
r  >  0.70 

RRI  HF 

.55  ±  .09 

44% 

RRI  LF 

.32  ±  .07 

0% 

RRI  HFnu 

.65  ±  .08 

55% 

RRI  LFnu 

.66  ±  .07 

55% 

RRISD 

.64  ±  .07 

44% 

pNN50 

.65  ±  .06 

55% 

Values  are  means  ±  SE.  MSNA,  muscle  sym¬ 
pathetic  nerve  activity;  RRI  HF,  R-R  interval 
high  frequency  power;  RRI  LF,  R-R  interval  low 
frequency  power;  RRI  HFnu,  normalized  high 
frequency  power,  RRI  LFnu,  normalized  low 
frequency  power;  RRISD,  R-R  interval  standard 
deviation;  pNN50,  the  percentage  of  adjacent 
R-waves  that  vary  by  50  msecs  or  more. 


we  chose  a  sample  cohort  of  subjects  who 
provided  MSNA  data  up  to  -70  mm  Hg 
(n  =  9)  and  ran  individual  correlations 
among  primary  associations  of  interest.  We 
chose  this  sample  cohort  because  including 
data  through  -70  mm  Hg  provided  a  min¬ 
imum  number  of  data  sets  for  robust  cal¬ 
culation  of  individual  correlation  coeffi¬ 
cients.  We  considered  a  correlation  of  >.7 
to  be  strong.  Average  individual  correlations 
with  the  percentage  of  individuals  with  cor¬ 
relations  of  >.7  are  shown  in  Table  2. 

RRISD,  pNN50,  RRIHF,  and  normalized 
RRIhf,  correlated  directly  with  decreases  of 
pulse  pressure  during  LBNP  (all  R2  >  .88). 
Normalized  RRILF  correlated  inversely  with 
pulse  pressure  (R2  =  .94),  but  absolute 
RRIlf  did  not  (R2  =  .005). 

DISCUSSION 

We  submitted  a  group  of  healthy  hu¬ 
man  volunteers  to  progressive  reduction 
in  central  blood  volume  to  determine 


whether  heart  period  variability  provides 
insight  into  autonomic  responses  to  pro¬ 
gressive  central  hypovolemia.  The  pri¬ 
mary  new  findings  from  this  study  are:  1) 
time-  (RRISD,  pNN50)  and  frequency- 
(RRIhf)  derived  components  of  heart  pe¬ 
riod  variability  decrease  progressively 
with  increases  of  MSNA,  and  2)  RRILF 
does  not  change.  Our  results  have  impli¬ 
cations  for  monitoring  of  bleeding  pa¬ 
tients  using  a  simple  ECG  and  provide 
insight  into  the  utility  of  time  and  fre¬ 
quency  domain  analyses  of  cardiac  inter¬ 
beat  intervals  as  markers  of  autonomic 
and  hemodynamic  compensation  to  cen¬ 
tral  blood  volume  reductions  and  early 
indicators  of  hemorrhage  severity. 

The  report  of  the  Post  Resuscitative 
and  Initial  Utility  of  Life  Saving  Efforts 
(PULSE)  trauma  work  group,  sponsored, 
in  part,  by  the  National  Heart,  Lung,  and 
Blood  Institute,  revealed  that  hemor¬ 
rhage  accounts  for  a  large  proportion  of 
deaths  due  to  traumatic  injury  (1).  Be- 
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cause  a  full  40%  of  traumatically  injured 
patients  die  before  they  reach  a  hospital, 
one  recommendation  of  the  work  group 
was  to  improve  or  develop  new  ap¬ 
proaches  to  monitoring  bleeding  patients 
(1).  Current  standard  vital  sign  monitors 
include  an  ECG,  periodic  blood  pressure 
recordings,  and  estimates  of  arterial  oxy¬ 
gen-hemoglobin  saturation.  The  need 
for  a  lifesaving  intervention  has  been  as¬ 
sociated  with  systolic  blood  pressure  of 
<90  mm  Hg  (24,  25);  unfortunately,  by 
the  time  these  pressure  reductions  occur, 
it  may  already  be  too  late  to  effectively 
resuscitate  a  bleeding  patient.  Cardiovas¬ 
cular  collapse  can  occur  quickly,  and 
without  advanced  warning  from  standard 
vital  signs.  Because  the  autonomic  ner¬ 
vous  system  is  fundamentally  responsible 
for  maintaining  stable  hemodynamics,  it 
is  appropriate  to  focus  on  autonomic  re¬ 
sponses  that  change  early.  The  results  of 
the  present  investigation  suggest  that 
RRISD,  pNN50,  and  RRIHF  may  provide 
newly  derived  vital  signs  that  could  be 
obtained  from  ECG  waveforms  available 
on  current  medical  monitors  to  assess 
the  need  for  early  intervention  in  advance 
of  alterations  in  standard  vital  signs. 

In  addition  to  hemorrhage,  it  is  likely 
that  injury  and  pain  associated  with 
trauma  may  result  in  highly  variable 
heart  rate  responses  to  trauma,  unrelated 
to  blood  loss  per  se,  hence  reducing  the 
sensitivity  of  heart  rate  as  a  predictor  of 
patient  outcome.  This  view  is  supported 
by  a  recent  analysis  of  10,825  patients 
with  blunt  or  penetrating  trauma,  reveal¬ 
ing  that  heart  rate  alone  was  neither  sen¬ 
sitive  nor  specific  in  determining  the 
need  for  emergent  intervention  (26).  Al¬ 
though  heart  rate  had  an  amalgamated 
correlation  of  0.83  with  MSNA  during 
progressive  LBNP  in  our  subjects,  RRIHf 
correlated  higher  (0.98)  with  MSNA  (Fig. 
3)  and  changed  earlier  (LBNP  at  -30  mm 
Hg)  (Fig.  2)  than  heart  rate  (LBNP  at 
-45  mm  Hg)  (Table  1)  during  central 
hypovolemia.  Our  results  therefore  sug¬ 
gest  that  an  index  of  sympathetic  nerve 
activity  may  be  more  sensitive  for  assess¬ 
ment  of  blood  loss  in  trauma  patients 
than  heart  rate.  This  hypothesis  is  sup¬ 
ported  by  the  observation  of  significant 
differences  in  heart  period  variability  be¬ 
tween  trauma  patients  who  died  com¬ 
pared  with  a  cohort  of  patients  who  lived, 
despite  their  heart  rates  being  similar  (7, 
8,  27).  Our  results  from  the  present  in¬ 
vestigation  are  therefore  consistent  with 
those  data  obtained  from  the  trauma  lit¬ 
erature  suggesting  that  the  use  of  heart 


period  variability  may  increase  the  sensi¬ 
tivity  of  decision-support  tools  for  clinical 
assessment  of  bleeding  patients. 

Assessed  in  the  time  domain,  increases 
or  decreases  of  heart  period  variability  are 
thought  to  reflect  increases  or  decreases 
of  the  overall  magnitude  of  autonomic 
influences  on  the  heart,  but  contribu¬ 
tions  from  either  parasympathetic  or 
sympathetic  components  cannot  be  de¬ 
termined  (22).  Nevertheless,  based  on 
thousands  of  individual  bedside  records 
from  patients  monitored  continuously  in 
the  ICU,  reductions  of  heart  rate  standard 
deviations  assessed  in  bins  of  1  to  4  secs 
are  associated  strongly  with  poor  out¬ 
come  (5,  6,  28).  Whether  poor  patient 
outcomes  assessed  with  periodic  analysis 
of  heart  rate  standard  deviations  are  as¬ 
sociated  specifically  with  autonomic  dys¬ 
function  has  not  been  determined.  In  this 
regard,  our  experiment  is  unique  in  that 
it  provided  the  opportunity  to  assess  the 
relationship  between  heart  period  vari¬ 
ability  derived  by  time  domain  and  direct 
measures  of  increasing  sympathetic 
nerve  activity  in  unanesthetized  humans 
undergoing  progressive  reduction  in  cen¬ 
tral  blood  volume.  Our  results  support 
the  notion  that  periodic  analysis  of  heart 
rate  standard  deviations  can  be  used  to 
assess  autonomic  dysfunction  in  bleeding 
patients  and  signal  the  need  for  early 
intervention. 

Although  autonomic  cardiovascular 
control  cannot  be  assessed  directly  in 
trauma  patients  (i.e.,  direct  measures  of 
MSNA),  the  results  of  our  experiment 
suggest  that  autonomic  function  may  be 
estimated  in  casualties  with  mild  to  se¬ 
vere  hemorrhage  through  analysis  of 
heart  period  variability  in  the  frequency 
domain.  In  humans,  cardiac  interbeat  in¬ 
tervals  oscillate  at  two  primary  frequen¬ 
cies;  high  (respiratory)  frequencies  (0.15- 
0.4  Hz)  and  low  frequencies  (0.04-0.15 
Hz).  High-frequency  oscillations  are 
eliminated  with  large  doses  of  atropine 
(29,  30),  and  low-frequency  oscillations 
are  reduced  with  propranolol  (31,  32). 
Based  largely  on  these  autonomic  block¬ 
ade  studies,  RRIHf  oscillations  have  been 
attributed  primarily  to  fluctuations  of  va¬ 
gal  nerve  traffic,  and  RRILF  oscillations 
have  been  attributed  to  a  combination  of 
factors,  including  fluctuations  of  both  va¬ 
gal  and  sympathetic  traffic  (9).  Although 
RRIlf  oscillations  are  associated  in  some 
way  with  fluctuations  of  sympathetic 
neural  activity  (10,  33),  it  is  clear  that  the 
association  is  imperfect  (14,  15).  Our  re¬ 
sults  support  this  concept  and  challenge 


the  usefulness  of  RRILF  as  a  tool  to  assess 
patient  status  during  hemorrhage. 

However,  even  those  investigators 
who  have  championed  the  use  of  spectral 
analysis  as  an  indicator  of  autonomic 
function  concede  that  RRILF  is  most  use¬ 
ful  as  an  estimate  of  sympathetic  nerve 
activity  during  procedures  that  progres¬ 
sively  increase  sympathetic  traffic  (9,  33, 
34).  Our  results  in  the  present  study  do 
not  support  this  notion  because  LBNP 
failed  to  increase  RRILF  in  even  a  single 
subject,  despite  progressive  increases  of 
MSNA  (Figs.  2  and  3;  0%  significant  cor¬ 
relations  between  RRILF  and  MSNA;  Table 
2).  It  therefore  seems  likely  that  sympa¬ 
thetic  compensation  to  central  hypovole¬ 
mia  as  occurs  during  hemorrhage  cannot 
be  adequately  estimated  or  tracked  di¬ 
rectly  with  RRIlf-  After  normalization  of 
RRIlf  to  express  LF  oscillations  as  a  frac¬ 
tion  of  total  power,  normalized  RRILF  in¬ 
creases  directly  with  MSNA  during  LBNP. 
However,  because  the  normalization  pro¬ 
cedure  involves  dividing  LF  by  LF+HF, 
and  we  observed  that  RRIHf  oscillations 
track  changes  of  directly  measured  MSNA 
whereas  RRILF  changes  do  not,  the  fact 
that  normalized  RRILF  tracks  MSNA  may 
simply  reflect  the  influence  of  the  rela¬ 
tionship  between  RRIHf  and  MSNA. 
Therefore,  our  experiment  may  be  the 
first  to  demonstrate  that  RRIHF  repre¬ 
sents  the  most  specific  frequency-derived 
heart  period  variability  measure  to  reflect 
the  autonomic  response  to  hemorrhage. 

Autonomic  responses  to  simulated 
hemorrhage  in  our  cohort  of  healthy  sub¬ 
jects  were  appropriate  as  evidenced  by  a 
relative  lack  of  change  in  mean  arterial 
pressure.  Compared  with  baseline,  sys¬ 
tolic  pressure  was  reduced  at  the  -45 
mm  Hg  level,  but  the  magnitude  of 
change  was  small  (about  7%),  and  sys¬ 
tolic  pressure  never  fell  to  <112  mm  Hg 
on  average,  despite  large  reductions  of 
central  volume  as  estimated  from  pulse 
pressure  (Table  1).  In  contrast,  as  early  as 
the  —30  mm  Hg  level,  both  pNN50  and 
RRIhf  were  reduced  by  about  50%  com¬ 
pared  with  baseline  (Fig.  2).  These  data 
highlight  the  fact  that  autonomic  func¬ 
tion  changes  earlier  than  standard  vital 
signs,  such  as  arterial  pressure  during 
simulated  hemorrhage,  and  suggest  that 
tracking  measures  of  autonomic  function 
may  provide  unique  insight  into  the  clin¬ 
ical  status  of  bleeding  patients.  Our  data 
further  indicate  that  the  application  of 
RRIhf  as  a  potential  clinical  tool  for  early 
assessment  of  blood  loss  compared  with 
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measurements  of  standard  vital  signs  can 
be  applied  universally  to  both  male  and 
female  patients. 

Without  substantial  anxiety,  pain,  or 
tissue  injury,  we  are  able  to  assess  with 
confidence  the  isolated  effect  of  progres¬ 
sive  central  hypovolemia.  Our  model  of 
simulated  hemorrhage  necessarily  in¬ 
duces  early  physiologic  compensations  to 
moderate  or  severe  central  hypovolemia, 
with  severe  hypotension  observable  at 
only  the  highest  levels  of  negative  pres¬ 
sures  and  just  before  cardiovascular  col¬ 
lapse.  In  some  severely  injured  patients, 
relative  hypotension  in  the  prehospital 
arena  is  sufficient  to  recommend  imme¬ 
diate  intervention  (24),  but  by  the  time 
arterial  pressures  fall  to  dangerous  levels, 
some  patients  may  be  unsalvageable.  The 
inability  to  maintain  adequate  arterial 
pressures  could  be  related  importantly  to 
autonomic  decompensation.  For  exam¬ 
ple,  Converse  et  al.  (35)  showed  that  hy¬ 
potension  induced  by  hemodialysis  was 
associated  with  a  relative  parasympa¬ 
thetic  activation  and  sympathetic  with¬ 
drawal.  Others  have  also  reported  para¬ 
sympathetic  predominance  (estimated 
from  an  elevated  HF/LF)  in  ICU  patients 
who  eventually  died  of  their  injuries  (3, 
4).  In  previous  studies,  we  estimated  au¬ 
tonomic  function  from  heart  period  vari¬ 
ability  in  prehospital  trauma  patients  and 
found  an  elevated  HF/LF  in  those  patients 
who  either  eventually  died  of  hemorrhage 
(7)  or  from  all-cause  trauma  (8).  In  all 
these  cases,  an  increase  in  parasympa¬ 
thetic  predominance  estimated  from  an 
elevation  in  RRIHf  was  associated  with 
increased  mortality  (i.e.,  a  late  indicator 
of  outcome).  The  results  from  our 
present  experiment  are  the  first  to  sug¬ 
gest  that  progressive  reductions  in  RRIHf 
can  provide  an  early  indication  of  the 
sympathetic  predominance  associated 
with  appropriate  compensation  during 
reduction  in  central  blood  volume  before 
reaching  a  condition  of  irreversible  pa¬ 
tient  outcome. 

Our  results  also  confirm  that  pulse 
pressure  may  provide  insight  into  actual 
volume  loss  because  it  tracked  the  reduc¬ 
tion  in  central  blood  volume.  Convertino 
et  al.  (21)  have  shown  that  arterial  pulse 
pressure  provides  a  good  indication  of 
relative  reductions  in  stroke  volume  dur¬ 
ing  central  hypovolemia.  More  impor¬ 
tant,  our  data  suggest  that  indices  of  au¬ 
tonomic  responses  (e.g.,  RRIHF,  RRISD, 
and  pNN50)  may  provide  noninvasive  de¬ 
cision  support  for  actual  circulating 
blood  volume  status.  Because  alterations 


in  RRIHF,  RRISD,  and  pNN50  are  highly 
correlated  with  pulse  pressure  (all  R 2  & 
0.88),  we  hypothesized  that  heart  period 
variability  metrics  should  represent  a 
sensitive  and  specific  marker  to  track  the 
effectiveness  of  fluid  resuscitation  and 
blood  loss.  Our  LBNP  model  provided  the 
opportunity  to  test  this  hypothesis  be¬ 
cause  central  blood  volume  was  restored 
to  baseline  at  the  cessation  of  negative 
pressure.  Indeed,  restoration  of  the  mag¬ 
nitude  of  heart  period  variability  (Fig.  2) 
with  concomitant  restoration  of  pulse 
pressure  (Table  1)  at  the  cessation  of 
LBNP  suggests  how  heart  period  variabil¬ 
ity  might  also  prove  useful  in  monitoring 
adequate  resuscitation  during  infusion  of 
fluids.  Thus,  in  addition  to  measurements 
of  blood  pressure,  continuous  monitor¬ 
ing  of  heart  period  variability  may  pro¬ 
vide  medical  personnel  insight  into  the 
status  of  central  blood  volume  during 
bleeding  and  intervention. 

Results  of  the  current  study  should  be 
interpreted  with  at  least  one  technical 
caveat,  and  that  is  the  practical  utility  of 
assessing  heart  period  variability  in  real 
time.  Due  to  the  necessity  of  identifying 
R-wave  fiducial  points  from  the  ECG  with 
absolute  certainty  for  frequency  domain 
analyses,  and  given  any  confounding  fac¬ 
tors  such  as  movement  artifact,  ectopic 
beats,  or  electrically  noisy  ECG  record¬ 
ings,  R-R  interval  frequency  domain  anal¬ 
yses  may  not  be  possible  for  clinical  ap¬ 
plications  until  future  monitors  have  the 
computing  capability  to  filter,  clean,  and 
distinguish  sinus  node  R-waves  from 
non-sinus  node  sources.  On  the  other 
hand,  it  seems  from  our  analyses  that 
simple  time  domain  heart  period  variabil¬ 
ity  metrics  also  track  autonomic  function 
well.  Such  analyses  are  more  forgiving  of 
ectopic  or  missing  beats  and  may  at  this 
time  be  the  most  promising  approach  to 
assess  autonomic  function  in  the  prehos¬ 
pital  setting  similar  to  what  has  been 
demonstrated  in  the  ICU  (5,  6,  28). 

Study  Limitations.  We  do  not  suggest 
that  subjects  in  the  present  study  resem¬ 
ble  actual  bleeding  patients  approaching 
hemorrhagic  shock.  Subjects  in  our  ex¬ 
periment  were  not  anxious  or  in  pain,  nor 
did  they  suffer  tissue  injury,  and  we 
therefore  cannot  speculate  how  the  addi¬ 
tion  of  these  real-world  issues  might  af¬ 
fect  our  results.  Although  we  are  unaware 
of  any  evidence  to  suggest  that  injury  and 
pain  would  alter  the  fundamental  physi¬ 
ologic  relationship  between  heart  period 
variability  and  sympathetic  nerve  activity 
observed  in  this  study,  we  cannot  dismiss 


the  possibility  that  trauma  in  addition  to 
central  blood  loss  could  influence  this 
relationship.  A  second  potential  limita¬ 
tion  of  the  current  study  is  the  loss  of 
nerve  recordings  at  high  levels  of  nega¬ 
tive  pressure.  Due  to  inherent  intricacies 
associated  with  obtaining  direct  measure¬ 
ments  of  MSNA  in  humans,  coupled  with 
movements  of  the  microelectrode  associ¬ 
ated  with  LBNP,  it  is  difficult  to  maintain 
adequate  nerve  recordings  as  chamber 
pressures  increase.  Despite  this,  to  our 
knowledge,  the  MSNA  data  we  present  at 
-70  and  -80  mm  Hg  are  the  only  data 
available  that  provide  insight  into  associ¬ 
ations  of  heart  period  variability  and  sym¬ 
pathetic  traffic  during  central  hypovole¬ 
mia  of  a  magnitude  similar  to  severe 
hemorrhage. 

Clinical  Implications.  Monitoring  and 
assessing  circulatory  volume  status  in 
bleeding  patients  from  standard  vital 
signs  is  problematic  because  vital  signs 
change  relatively  late  during  injury  pro¬ 
gression.  Monitoring  of  heart  period  vari¬ 
ability  may  increase  the  capability  to 
identify  patients  who  are  or  are  not  com¬ 
pensating  appropriately  to  hemorrhage. 
Potential  uses  include  remote  monitor¬ 
ing  in  mass-casualty  situations,  prehospi¬ 
tal  emergency  transport  (7,  8),  ICU  mon¬ 
itoring  of  septic  patients  (36)  and  patients 
undergoing  renal  dialysis  (37)  leading  po¬ 
tentially  to  hypotension,  and  other  con¬ 
ditions  in  which  blood  volume  has  the 
potential  to  decrease  undetected.  In  the 
prehospital  or  military  setting,  tracking 
changes  of  autonomic  function  with 
heart  period  variability  from  a  simple 
ECG  could  increase  the  potential  to  rec¬ 
ognize  a  bleeding  patient  who  seems  sta¬ 
ble  but  who  may  actually  be  progressing 
to  cardiovascular  collapse  and  death. 
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